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. INTRODUCTION 

O 

^ ■ It is currently believed that the gravitational waves that come from the final 

(3J[). stages of the evolution of compact binaries just before their coalescence are very 

likely signals to be detected by long- arm laser interferometers M. The two projects 
to build such detectors - LIGO and VIRGO - are already approved and they are 
rapidly progressing, and the third one - GEO600 - is likely to be funded soon. A 
standard optimal method to detect the signal from a coalescing binary in a noisy data 
set and to estimate its parameters is to correlate the data with the filter matched 
to the signal and vary the parameters of the filter until the correlation is maximal. 
It has recently been realized^ that the correlation is very sensitive even to very 
small variations of the phase of the filter because of the large number of cycles in the 
signal. Consequently, the addition of small corrections to the phase of the signal due 
to the post-Newtonian effects decreases the correlation considerably. Thus the post- 
Newtonian effects in the coalescing binary waveform canbe detected and estimated 
to a much higher accuracy than it was thought beforeH. We present an analysis 
of the estimation of parameters of the post-Newtonian signal. We also examine the 
detectability of the post-Newtonian signal and estimation of its parameters using 
the Newtonian waveform as a filter. This filter can be used as the simplest search 
template. 

GRAVITATIONAL- WAVE SIGNAL FROM A BINARY 

Let us first give the formula for the gravitational waveform of a binary with 
the currently known post-Newtonian corrections. In this paper we work within 
the so called "restricted" post-Newtonian approximation, i.e., we only include the 

1 This work was supported by KBN Grant No. 2 P302 076 04 



post-Newtonian corrections to the phase of the signal keeping the amplitude in its 
Newtonian form; this is because the effect of the phase on the correlation is domi- 
nant. The inclusion of post-Newtonian effects in amplitudes will not qualitatively 
change our results. Due to radiation reaction the orbit of the binary is rapidly cir- 
cularized, nevertheless, we include the first order correction due to eccentricity for 
completeness^. The tidal effects are known to be very small and we neglect them. 
We also include the contributions due to dipole radiation predicted by the Jordan- 
Fierz-Brans-Dicke (JFBD) theory (assuming that they are small) to investigate the 
possibility of testing GR against alternative theories of gravity. 

The analysis of the signal is best performed in the Fourier domain. The expres- 
sion for the Fourier transform of our signal in the stationary phase approximation 
is given by 
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where e is the eccentricity of the binary at gravitational frequency f (k e is constant 
to order e 2 ,), s Q and s s are spin-orbit and spin-spin parameters, respectively^ . The 
constants c% and c 2 in the 2nd post-Newtonian contribution are only recently being 
calculated || and calculation of further corrections is in progress. Cj,i=(l,2), is 
the sensitivity of the body i to changes of the scalar field, F(Ci, C 2 , mi, m 2 ) is a 
complicated function of sensitivities and masses of the bodies and to is the parameter 
of the JFBD theory. Current observational tests restrict u > 600 and observations 
of the Hulse- Taylor binary pulsar restrict u> > 200. 

ESTIMATION OF THE POST-NEWTONIAN PARAMETERS 



A basic method proposed to detect the gravitational-wave signal of a coalescing 
binary and to estimate its parameters is matched-filtering and maximum likelihood 

2 This correction was derived by N. Wex 



(ML) estimation. The performance of the method is determined by signal-to-noise 
ratio d and the Fisher information matrix T. Signal-to- noise ratio determines the 
probability of detection of the signal and the diagonal elements of the inverse of the 
Fisher matrix are approximately the variances of the maximum likelihood estimators 
of the parameters for large d. 

We have considered three representative binary systems with parameters sum- 
marized in Table I. 

Table I. Parameters of Neutron Star (NS) and Black-Hole (BH) Binary Systems 



Binary 


miM 


m 2 M 


MM Q 


s 


kM~ 5/3 




h/2M & 2/3 


NS-NS 


1.4 


1.4 


1.2 


2.4xl0" 2 


0.72 


4.1 


25 


NS-BH 


1.4 


10 


3.0 


4.0 


0.16 


2.0 


16 


BH-BH 


10 


10 


8.7 


3.9 


2.7xl0" 2 


0.58 


4.7 



Assuming that the noise is Gaussian, neglecting the effects of eccentricity, 2nd 
post-Newtonian effects, dipole radiation, and assumining that binaries are located at 
a distance of 200Mpc we have calculated the signal-to-noise ratios and approximate 
values of the variances of the ML estimators. The results are summarized in Table 
II (cf. § Table II). 

Table II. Signal-to-Noise Ratio and RMS Errors for Parameters of the 3 Binaries 



Binary 


S/N 


s/n 


At a ms 


AM/M 




Am/m 


As/s 


NS-NS 


15 


32 


0.75 


0.023% 


6.4% 


9.6% 


33 x W 2 % 


NS-BH 


32 


15 


0.65 


0.061% 


4.8% 


7.0% 


8.2% 


BH-BH 


77 


6 


0.46 


0.19% 


16% 


23% 


38% 



The number y/n gives the improvement of the signal-to-noise ratio due to matched 
filtering. 

We have also examined the effects of eccentricity and dipole radiation We have 
added their contributions to the phase of the 3/2post-Newtonian signal separately 
and we have calculated the Fisher matrix and its inverse. For a NS-NS binary the 
relative rms error in eccentricity is given by — = 0,63x10 — . For the Hulse - Taylor 



binary pulsar we have a e = e^f 



2 /19/9 



1.8 x 10 13 . Thus the effects of eccentricity 



would be practically undetectable for such a binary. For the NS-BH binary the 
relative rms error in the dipole radiation coefficient ku in JFBD theory is given by 
0.75 (ygo)- Thus, the effects of the dipole radiation could be determined to 
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about the same accuracy as from timing of the binary pulsar in our Galaxy. 

It is well known that the rms errors of the estimators of the parameters increase 
with the number of parameters. We have investigated this effect with the increasing 
number of post-Newtonian parameters. We have considered a reference binary of 
M = IMq located at the distance of lOOMpc. 

Table III RMS Errors for Parameters at Various post-Newtonian Orders 



At c msec 


A<f> c 


AkM~J>' 3 


AhM' 1 


Ak 3/2 M~ 2/3 


Ak 2 M~ 1/3 


Ak e M~ w/9 


0.17 


0.10 


8.3 xl0~ e 










0.27 


0.33 


4.0 xlO" 5 


5.8xl0~ 3 








0.54 


1.9 


1.7 xlO" 4 


0.70 xlO- 1 


0.52 






1.6 


24 


6.6 xl0~ 4 


0.50 


7.2 


28 




2.3 


45 


2.3 xl0~ 3 


1.3 


17 


59 


1.2xl0" 6 



We were able to include the 2nd post-Newtonian correction because the Fourier 
transform of the signal is linear in the mass parameters ki and consequently, the 
Fisher matrix for these parameters does not depend on their numerical values. In 
Table IV we show the degradation of accuracy of estimation of the chirp mass, the 
reduced mass, and the total mass with the increasing number of parameters in the 
template for the NS-NS binary at a distance of 200Mpc. 

Table IV. RMS Errors for Masses at Various post-Newtonian Orders 



pN order 


AM/M 


A/i//x 


Am/m 


1 pN 
3/2 pN 

2 pN 


0.0054% 
0.023% 
0.080% 


0.55% 
6.4% 
42% 


0.81% 
9.6% 
63% 



The rms errors in [i and m do not depend on the value of 2nd post-Newtonian 
mass paramter ki- 

THE NEWTONIAN FILTER 

On the one hand the correlation of the signal with the template is very sensitive 
to small corrections in the phase of the signal on the other hand the accuracy of 
estimation of the parameters is significantly degraded with increasing number of 
corrections even though a correction may be small. Moreover we cannot entirely ex- 
clude unpredicted small effects in the gravitational-wave emmission (e.g. corrections 
to general theory of gravity) that we present cannot model. Thus, there is a need 
for simple filters or search templates that will unable to scan the data effectively and 
isolate stretches of data where the signal is most likely to beo . The simplest such 
filter is just a Newtonian waveform hw the Fourier transform of which, in stationary 
phase approximation, is given by Iin = / _7 / 6 expi[27r/i c — 4> c — 7r/4 + fc^|g(7r/) -5 / 3 ]. 
Using such a suboptimal filter decreases the correlation and, consequently, decreases 
the signal-to-noise ratio. Also the estimate of the k parameter is shifted by a definite 
amount depending on the noise of the detector and the parameters of the binary. 
In Table V we have given the drop in signal-to-noise given by FF = (h\h^)/(h\h) 
where (/11I/12) is the correlation of functions hi and /12 (see Ref.|| Eq.(2.3) for defini- 
tion) and the parameter I = y/~FF. In the case of the Newtonian filter probability of 
detection is determined by two parameters: optimal signal-to-noise ratio d = \J (h\h) 
and d a = a/ (h\h ) = I x d. We have also calculated the shift in the estimate of the 
k parameter and the accuracy of the determination of k. We have included the 1st 
and the 3/2 post-Newtonian corrections. 



Table V. Performance of the Newtonian Filter 



Binary 


I 


FF 


Shift 5kM~ 5/3 


Accuracy A/cM Q 5 ^ 3 


NS-NS 


0.90 


0.81 


0.01560 


0.13 x 10" 3 


NS-BH 


0.87 


0.76 


0.004898 


0.026 x 10~ 3 


BH-BH 


0.98 


0.96 


0.001209 


0.011 x 10" 3 



When the amplitude and phase modulations are taken into account it was 
shown0 that in the worst case FF = 0.39 (1 = 0.63). 

Using the Newtonian filter we would not like to loose any signals. We can 
achieve this by suitably lowering the detection threshold when filtering the data 



with the Newtonian filter. By this procedure we would isolate stretches of data 
with all the signals that would be detected with optimal filter and also an increased 
number of false alarms. The next step would be to analyse the reduced set of data 
with accurate templates and the initial threshold to make the final detection and 
estimate the parameters of the signal. 

In Table VI we have given examples of the performance of the above procedure. 
We assume the detection threshold T = 5 and we assume that we have one signal 
for the optimal signal-to-noise ratio d. N is the expected number of detected signals 
with the optimal filter, Np is number of false alarms, Nn is the number of detected 
signals with the Newtonian filter, T\r is the lowered threshold, Nl is number of 
signals with the lowered threshold and Nfl is number of false alarms with the 
lowered threshold. 



Tab 


le VI 


d 


FF 


N 


N F 


N N 


T N 


N L 


N FL 


15 


.81 


27 


0.055 


20 


4.5 


28 


0.16 


15 


.36 


27 


0.055 


5.6 


3.225 


28 


2.1 


30 


.81 


225 


1.1 


165 


4.5 


230 


2.2 


30 


.25 


225 


1.1 


31 


2.875 


229 


32 



A different search template consisting of the post-Newtonian waveform with spin 
effects stripped off has been analysed in Ref . J5J . 
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